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DISCLAIMER 
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warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
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DESCRIPTION: The solution to the initial value problem 



- f (x,y), y(x Q ) - y Q 

is numerically approximated and solution pairs (x,y) are printed 
on the typewriter at a specified printing increment, HPR, from 
x=x q to x-XEND. The solution can also be plotted on the 

typewriter if desired. MOTE: x < XEMD. 

The Runge-Kutta-Gill Method may be used (HRK » 0.0) or the 
Runge-Kutta-Gill Method may be used to furnish starting values 
for Hammings Method (KRK ■ 1.0). For a description of these 
methods rf»f pr to Mathematical Methods for Digital Computers 
by A. Ralston and H. Wilf . * 

The stepping increment is denoted by H and the solution value 
by Yl. 

When computing by Hammings Method, three numbers, Tl, T2, and 
T3» must be provided by the user. These numbers control 
truncation error as described in the listing. 



INPUT/OUTPUT: Input may be by cards or typewriter. Sample input for cards 
is included i n the listing. 

Output is on the typewriter and four pairs of sample outputs 
are included (one for Runge-Kutta-Gill and one for Hammings). 



SAMPLE PROBLEM: Four sample problems are described in the listing. In order to 
use the program deck, statement 100, which contains no GOTRAN 
statement, must be pulled out and replaced by the relevant 
function computation^ samples of which are given in the listing. 

Operating instructions are given in the listing. 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 




PROGRAM LISTING 



C GOTRAN WITH .SUBROUTINE TECHNIQUE USFD FOR HAMMINGS METHOD AND /OR 
C RUNGE KUTTA GILL METHOD FOR INITIAL VALUE PROB INVOLVING ONE 
C FIRST ORDER D IFF . FO. 

C RFFFR TO GOTRAN OPFRATING MANUAL FOR TYPEWRITER TAB SETTINGS 
C PROGRAM STOPS ON RFTURN TO ^TATFMFNT 3 
C HRK=0.0. RKG ONLY — , =1.0, HAM WITH RtCG 

C IN HAM. T LESS THAN T2 DOUBLE H, T GREATFR THAN Tl HALVE H, T LESS 
C THAN T3 PRINT O.n AND RFTURN TO 3 AND STOP 
C PLT=O.P, NO PLOT, =1.0, SOLUTION IS 

C TRANSLATED BY YAXS AND THIS RESULT IS SCALED BY YSC AND THIS LAST 
C RESULT PLOTTED ON THE TYPEWRITER. CHOOSE YSC SO t Yl +YAXS >/YSC IS IN 
C TO 49. IF NOT POSSIBLE, SOLUTION WILL BE PRINTED BUT NOT PLOTTED 
C WHEN (Y1+YAXS)/YSC IS NOT IN THIS RANGE. 

DIMENSION YH{3> 
3 READ ,HRK,DLT 

R C AD ,HPR,XFND,H 

PRINT ,H 
XPR=HPR 
RF An ,X0,Y10 

IF (HRK) 31,30,31 
31 READ ,T1,T2,T3 

PRINT ,T1,T2,T3 
30 PRINT ,X0,Y10 

C°=0.0 

IF (OLTK301 ,301 ,900 
C END DATA INPUT 
C BEGIN RKG 
301 X = XO 

Y1=Y10 

01 =0.0 

K=l 

CK=1 .0 

I U = X 

V = Y1 
C I=! . 

GO TO 100 

II C1=H*F1 
TEMP=2»0*Q1 
TEMP=C1-TEMP 
Dl = .°i»TEMP 
Wl =Y1 +D1 
TFMP=.5*C1 
TFM1 s3,*01 
TFMP = TP"M 1 -TEMP 
Ol =01 +TEMP 
Us.^+H 

U = X+U 

V = W1 
CI=2. 

GO TO 100 
18 C1=H*F1 

TFMP=C!-01 

01 = .?9?BQ3? c ;*TFMP 

Wl =W1 +D1 

TFMPr • 2P2B932 C ;*C 1 
TFM 1 =7.0*01 



TFMP.TEM1 -TEMP 
01«01+TEMP 
V»W1 
CI*3.0 
GO TO 100 
25 C1«H#FI 

01 ■! .7071 067#TEMP 

wi«wi+ni 

TFMP-1 •7071067*C1 
T*Mt ■3,0*01 
TFMP-TEM1-TEMP 
OJ-OI+TCWP 
U"X+H 
V»W1 
C!«4.0 
60 TO 100 
32 C!»H#F! 

TFWP«2.0#Ql 

TF MP-C1-TEMP 

1 ■ • 1 6666667*TE WP 

TFWP-.5#C1 
TFM1«3«0*D1 
TE**P«TEM1 —TEMP 
01-Qt+TEMP 
X«X+M 

PPNT»X-XPP. 
IF (PWNT»2»200«200 
2 IF (HPK)lt!i3? 

37 YH(K) ■ VJ 

DK«CK-3.0 
TF (niO40«51 ( 5I 

CK«CK+I .0 

GO TO 1 
C END PKG 
C BEG IN HAMMING 
510 7I-0.0 

CK»4»0 

HJ«0.0 

Xl-XO+H 

Yt1«0.0 

Yl ?»0»0 

U-Xl 
J"l 

V » YH(J> 
CI-5.0 

GO TO 100 
513 F11«F1 
O-Xl+M 
J«2 

V ■ YH(J> 
C !»6.0 

GO TO 100 
Ff?»Fl 



Us?.0*H 

u»xi +u 

J*3 

V « YH(J) 
C 1*7.0 
GO TO 100 
522 F13»F1 

4 TEMP«2.0*FI3 
T FMP« TEMP-F 1 2 
TFM1 »2.0*F1 1 
TFMP-TEMP+TFM1 
TFMla! ,3333333*H 
TFMP«TFMP#TEM1 
PI «Y1 O+TEWP 

01 «.9?56 1983*21 
Q! *P1-Q1 
U«3.0#H 
U-Xl+U 
V«Q1 
Ct»8.0 
GO TO 100 
531 RlsFl 

TFMP»2.0*F13 
TFMP*TFMO+Pi 
TFWP«F1?-TFMP 
TFM1«3.0*H 
TFMPstTFM 1 *TEMP 
J«3 

TEM1 »9.0*YH(J) 
J*l 

TEM1 ■ YH(J)-TEM1 
TFMP«TEM1 +TEMP 
TEMP«.125*TEMP 
C1»0,0-TEMP 
71»P1-C1 

J»3 

TF (PI ) 534«535«535 

534 TY « 0.0 —PI 
GO TO 536 

535 TY »P1 
536 TEST-TY-1.0 

IF (TEST)537,537.S39 
537 T « .92561983*21 

GO TO 540 
539 T «01/P1 

T ■ .92561 983*T 
540 IF (T > 532 « 533 » 533 
53? T»0.0-T 
533 TFSTaT— T2 

TF (TEST J552.5.5 

5 TEST=T-T1 

IF (TEST ) 538*538*548 
538 Y13«Y12 
Yl 2*Yl 1 
Yl 1 = Y1 O 
J»l 




L=2 

M = 3 

Y10«YH(J) 

YH( J)«VH(L ) 

YH(L>«YH(M) 

YH(M>» .074380 165*Z1 

YH(M)«C1+YH(M) 
Fl 1*F12 
Fl 2«F!3 
X1«X1+H 
X»2»0#H 
X-Xl+X 
PPNT«X-XPR 
IF (PRNT>543«541 ,541 
•541 J«3 

Y1«YH( J) 
GO TO 200 
543 U«2»0»H 
U«X1+U 

J=3 

V*YH(J) 
CI-9.0 
GO TO 100 
546 F13=F1 
GO TO 4 

548 H«»5#H 

PRINT ,H 
TEST»H-T3 

IF (TEST 1549,550,550 

549 HM1N«0»0 
PRINT »HM IN 
GO TO 3 

•550 X0«4.0*H 
X0»X1+X0 
J=3 

Y10«YH( J> 
GO TO 301 
55? IF (HJ 1553,554. 553 

553 MJ«O.0 
GO TO 5 

554 TEST«4«0*H 

TEMP » XO + TEST 

TEM1 "TEMP+TEST 

TEST »X1-TEMP 

IF (TEST) 5*556.536 
556 TEST -TEM1-XEND 

IF (TEST ) 551,5,5 
551 H «2.0#H 

PRINT »M 
HJ«1 .0 
X0«2»0#H 
XOaXl -XO 

J"l 

L = ? 

YH(L)«YH(J) 
YH( J)»Y1 1 



s, 



Y10=Y13 

GO TO 510 
C END HAMMING 
C BEGIN PRINT 
200 PRINT ,X,Y1 

IF (PLT )201 .201 ,901 
?0 1 TEST * X-XENO 

IF (TFST)202,3,3 
202 XPR«XPR+HPR 

IF (HRK11 ,1,205 
205 TEST«CK-4»0 

IF (TEST)2,543,2 
C END PRINT 

C SUBROUTINE FOR PLOT — REMOVE IF COMPILED IN FORTRAN WITHOUT FORMAT 
C AND REPLACE BY 900 STOP AND 901 STOP 
900 READ »YAXS,YSC 

Y1«Y10 v 
pnj YP»Y1+YAXS 
YP»YP/YSC 
IF (YP1904.902 .902 

902 YT=YP-49» 

IF (YT)903,903,904 

903 PLOT (YP» . ) 

04 IF (CP1201 ,905,201 

05 CP*1.0 
GO TO 301 

C END PLOT 

C SUBROUTINE FOR F1*F(U»V> 

1 00 

zi»et-2.o 

IF (71 >i i .in.iot 

101 7!=C!-4,0 

IF (ZI 125,32, 1 02 

102 ZI=Ct-6.0 

IF (ZI 1513,518,103 

103 ZI=CI-8» 

IF (ZI 1522,531 ,546 
C END SUBROUTINE 
END 



SAMPLE PROBLEM DATA 



C THERE FOLLOWS 4 SAMPLE FUNCTION FORMULAS WHICH MUST RFPLACE THE 
C BLANK STATEMENT FOLLOWING 100 IN THE SOURCE DECK 

100 Fl« 1.0 

100 Fl» U 

100 Fl« V 
100 TFMP»t;»U7 
Fl «V#V? 
F1«F1+TEMPZ 
C SAMPLE DATA COMMON TO ALL FUNCTIONS USED 
C PROB. NO. I — RKG NO PLOT. 
0.0 
0.0 
0.1 
1 .0 

• 01 
0.0 
1 .0 

C PROB. NO. 2 — RKG WITH PLOT 

0.0 

1.0 

0.1 

1 .0 

• 01 
0.0 
1*0 
-1 .0 

• 0? 

C PROB. NO. 3 — HA MM IMG NO PLOT 

1 .0 

0.0 

0.1 

1 .0 

• 0!?«t 
0.0 

1 .0 

.00001 
.0000001 

• 001 

C PROB. NO. 4 — HAMMING WITH PLOT 
1 .0 
1 .0 

o.t 
i .o 

.01??! 

0.0 

1 .0 

.0001 

.0000001 

.001 

-1 .0 

.03 



I.OOOOOOOE- 2 
0.0 

.10000000 
.20000000 
.30000000 
.if 0000000 
.50000000 
.60000000 
.70000000 
.80000000 
.90000000 
1.0000000 
Hx I.OOOOOOOE- 2 
0.0 

.10000000 
.20000000 
.30000000 
.40000000 
.50000000 
.60000000 
.70000000 
.80000000 
.90000000 
1; 0000000 
H* 1.2500000E- 2 
1.0000000E- 5'Tl 
0.0 



1.0000000 
1.0999990 
1 . 1 999980 
1.2999970 
1.3999960 
1 .4999950 
1.5999940 
1 .6999930 
1 .7999920 
1.8999910 
1.9999900 

1.0000000 
1.0999990 
1.1999980 
1.2999970 
1.3999960 
U4999950 
1 ^5999940 
1.6999930 
P. 7999920 
1.8999910 
1.9999900 — *• "** 

1.0000000E- 7 *T g 
1.0000000 



i.oooooooe- 3 * T * 



2.5000000E- 


2 




.10000000 




1.0999997 


.20000000 




1.2000000 


.30000000 




1.3000005 


.40000000 




1.4000006 


*=5.0000000E- 


2 




.50000000 




1 .5000006 


#s .10300000 






.65000000 




1.6500006 


.75000000 




1.7500006 


.85000000 




1.8500006 


.95000000 




1.9500006 


1.0500000 




2.0500006 


0* 1.2500000E- 


2 




1.0000000E- 


4--TJ 


I.OOOOOOOE- 7 


0.0 




1.0000000 


fr«2.5000000E- 


2 




.10000000 




1.0999997 


.20000000 




1.2000000 


.30000000 




1.3000005 


.40000000 




1 .4000006 


H- 5.0000000E- 


2 




.50000000 




1.5000006 


ft- .10000000 






.65000000 




1.6500006 


.75000000 




1.7500006 


.85000000 




1 .8500006 


.95000000 




1.9500006 


1.0500000 




2.0500006 — * 



I.OOOOOOOE- 3 - 7> 



«*1.0000000E- 2 
0.0 

.10000000 

.20000000 

.30000000 

.40000000 

.50000000 

.60000000 

.70000000 

.80000000 

.90000000 

1.0000000 

1.0000000E- 2 

0.0 

.10000000 

.20000000 

.30000000 

* 4 0000000 

.50000000 

.60000000 

.70000000 

.80000000 

.90000000 

1.0000000 

1.2500OOOE- 2 

1.0000000E- 5»™ 

0.0 

H=6.2500000E- 3 

.10000000 
^ e 1.25 00000 E- 2 
2.5000000E- 2 



.21250000 




5. 0000000 E- 


2 


; 31 250000 




H*-. 10000000 




'.46250000 




.56250000 




.66250000 




W= .20000000 




.86250000 




1.0625000 




1.2500000E- 


2 


1.0000000E- 




0.0 




.10000000 




*s 2.5000000E- 


2 


.20000000 




«* 5. 000 0000 E- 


2 


.30000000 




.10000000 




.45000000 




.55000000 




.65000000 




.75000000 




.85000000 
.95000000 






.20000000 




1.1500000 





1.0000000 
1.0049160 
1.0198320 
1.0447480 
1.0796640 
1.1245800 
1.1794960 
1.2444120 
1.3193280 
1 .4042440 
1.4991600 

1.0000000 
1.0049160 
1.0198320 
1.0447480 
1.0796640 
1 . 1 245800 
1.1794960 
1.2444120 
1.3193280 
1.4042440 
1.4991600 

1.0000000E- 1'Tz 
1.0000000 

1.0053415 



1.0229196 

1.0491696 

1.1072946 
1.1585446 
1.2197946 

1.3722946 
1 .5647946 

1.0000000E- 7* T * 

1.0000000 

1.0049617 

1.0199616 

1.0449615 

1.1012115 
1.1512115 
1.2112116 
1.2812116 
1.3612116 
1.4512116 

1.6612116 




1.0000000E 



1.0000000E 



*» 1.0000000E- 2 
0.0 

. 1 0000000 
20000000 
.30000000 
.40000000 
.50000000 
'.60000000 
.70000000 
.80000000 
.90000000 

1; ooooooo 

fM.OOOOOOOE- 2 
0.0 

.10000000 
.20000000 
.30000000 
.40000000 
.50000000 
.60000000 
.70000000 
.80000000 
.90000000 
1 ".0000000 
U~ U2500000E- 2 
1.0000000E- &TL 

o-.o 

« = 6.2500000E- 3 

1.2500000E- 2 

; 10625000 

".20625000 

.30625000 
*» 2.5000000E- 2 

-.40625000 

.50625000 

5.0000000E- 2 

-.63125000 

-.73125000 

; iooooooo 

-.88125000 
.98125000 
1; 081 2500 
1.2500000E- 2 
T.OOOOOOOE- 4 = 72 
0.0 

.10000000 
. 20000000 
2.50O0OOOE- 2 
.31250000 
.41250000 
.51250000 
5.0000000E- 2 
.61250000 
.71250000 
.81250000 
.91250000 
1 .0125000 



1.0000000 
1.1050784 
1.2211980 
1 -.34951 91 
1.4913240 
1; 6480302 
1.8212027 
2.01 25722 
2.2240500 
2".4577497 
2.7160068 

1.0000000 

1.1050784 

1.2211980 

1.3495191 

1.4913240 

1.6480302 

1 -.821 2027 

2.0125722^ 

2; 2240500 

2;4577497 

2.7160068J 



UOOOOOOOE- l-.Tz 
1.0000000 



1; 11 20453 
1". 2290009 
1.3582566 

1 -.501 1058 
1.6589788 

1 - .8798694 
2.0775770 

2;41 37999 

2 - .66766 15 
2.9482220 

1.0000000E- 7-U 

r. ooooooo 

1.1051288 
1.2213570 

1.3667870 
1.5105335 
1.6693980 



1.8449701 
2.0390074^ 
2; 25345 18 
2-.4904495 V 
2.7523726 / 



4y s 



Co) 



1 .OOOOOOOE- 3 a TJ 



1.0000000E- 



c % r.oooooooE- 
o.o 

-.10000000 

-.20000000 
-.30000000 
-.40000000 

-.50000000 

'.60000000 
.70000000 
.80000000 
.90000000 
ERROR E3 
1; 0000000 
**1;0000000E- 
0.0 

.10000000 
'.20000000 
.30000000 
.40000000 

: 50000000 

'.60000000 

; 70000000 

'.80000000 
.90000000 

ERROR E3 

r. ooooooo 

tf* 1.250OOOOE- 
UOOOOOOOE- 
0;0 

6;2500000E- 
V* 3-.1250000E- 
"» 6.2500000E- 

.10000000 
"* 1.2500000E- 
'.20625000 
.30625000 
2.5000000E- 
.41875000 
1.2500000E- 
.50625000 
**6;2500000E- 
3.1250000E- 
H' 1-.5625000E- 
/y»3;1250000E- 
^6.2500000E~ 

.60468750 
H? 1.2500000E- 
.70468750 
6 .2500000E- 
*s 3.1250OO0E- 
I ^SOOOE- 
W.SUSOOOE- 

o-.o 

1.2500000E- 
KOOOOOOOE- 
0.0 

.10000000 
V* 2.5000000E- 
.21250000 
.31250000 
.41250000 
'.51250000 
.61250000 



2 

5«7v 

3 
3 
3 



1; OOOOOOO 
1.1112359 
1 -.2523885 
1 -.4383257 
1 -.693431 2 
2; 061 6388 
2.6325588 
3; 6246396 
5.7528798 
13.494457 

1. 3345425 E+42 



1; OOOOOOO 
1.11 12359 
1; 2523885 
1 '.4383257 
1 -.693431 2 
2; 061 6388 
2;6325588 
3*.6246396 
5.7528798 
13.494457 



1 .3345425 E+42 



4k 




1-.0000000E- 7*7-* 1.0000000E- 3 
1.0000000 



1.1113295 

i; 2629955 
1.4530976 

1.7547470 

2.0949455 



2.6764518 
3.7143019 



4»7i 1.0000000E- 7 *T* 1.0000000E- 3 * T3 

r. ooooooo 

1.1113658 

2 

1 .2733678 
1 -.4670806 
1 -.7347283 

2:7385827^ fit**^ ><- '«? 



10 



1.2500000E- 2 

.70000000 
H- 6:2500000E- 3 
M* 3.1250OOOE- 3 
^OSOOOOOE- 3 

-.80000000 

.90000000 

3-.1250000E- 3 

U5625000E- 3 
W i 7;8l25000E- 4 

0.0 



3*6506861 



5.8326228 
14.211229 



